December 24, 2009 13.32 



WSPC - Proceedings Trim Size: 9in x 6in ST-0083 



1 



DRIFT EFFECTS AND THE AVERAGE FEATURES OF 
COSMIC RAY DENSITY GRADIENT IN CIRS DURING 
SUCCESSIVE TWO SOLAR MINIMUM PERIODS 

A. FUSHISHITA 1 , Y. OKAZAKI 2 , T. NARUMI 1 , C. KATO 1 , S. YASUE 1 , T. 

KUWABARA 3 , J. W. BIEBER 3 , P. EVENSON 3 , M. R. DA SILVA 4 , A. DAL LAGO 4 , 

N. J. SCHUCH 5 , M. TOKUMARU 6 , M. L. DULDIG 7 , J. E. HUMBLE 8 , I. SABBAH 9 , 

J. KOTA 10 , and K. MUNAKATA 1 

1 Physics Department, Shinshu University, Matsumoto, Nagano 390-8621, Japan, 
2 Department of Geophysics, Tohoku University, Sendai, Miyagi 980-0861, Japan 
3 Bartol Research Institute and Department of Physics and Astronomy, University of 

Delaware, Newark, DE 19716, USA 
^National Institute for Space Research (INPE), 12227-010 Sao Jose dos Campos, SP, 

Brazil 

^Southern Regional Space Research Center (CRS/INPE), P.O. Box 5021, 97110-970, 

Santa Maria, RS, Brazil 
6 Solar Terrestrial Environment Laboratory, Nagoya University, Nagoya, Aichi 

464-8601, Japan 

7 Australian Antarctic Division, Kingston, Tasmania 7050, Australia 
8 School of Mathematics and Physics, University of Tasmania, Hobart, Tasmania 

7001, Australia 

9 Department of Physics, Faculty of Science, Kuwait University, Kuwait; on leave from 
Physics Department, Faculty of Science, University of Alexandria, Egypt 
10 Lunar and Planetary Laboratory, University of Arizona, Tucson, AZ 87721, USA 



We deduce on hourly basis the spatial gradient of the cosmic ray density in three 
dimensions from the directional anisotropy of high-energy (~50 GeV) galactic 
cosmic ray (GCR) intensity observed with a global network of muon detectors 
on the Earth's surface. By analyzing the average features of the gradient in 
the corotational interaction regions (CIRs) recorded in successive two solar 
activity minimum periods, we find that the observed latitudinal gradient (G z ) 
changes its sign from negative to positive on the Earth's heliospheric current 
sheet (HCS) crossing from the northern to the southern hemisphere in A<0 
epoch, while it changes from positive to negative in A>0 epoch. This is in 
accordance with the drift prediction. We also find a negative enhancement in 
G x after the HCS crossing in both A<0 and A>0 epochs, but not in G y . This 
asymmetrical feature of G x and G y indicates significant contributions from the 
parallel and perpendicular diffusions to the the gradient in CfRs in addition 
to the contribution from the drift effect. 
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1. Introduction 

The drift model of cosmic ray transport in the heliosphere predicts a bidirec- 
tional latitude gradient of the GCR density, pointing in opposite directions 
on opposite sides of the HCS [ x ][ 2 ]. The predicted spatial distribution of 
the GCR density has a minimum along the HCS in the "positive" polarity 
period of the solar polar magnetic field (also referred as the A>0 epoch), 
when the interplanetary magnetic field (IMF) directs away from (toward) 
the Sun in the northern (southern) hemisphere, while the distribution has 
the local maximum on the HCS in the "negative" period (A<0) with the 
opposite field orientation in each hemisphere. The field orientation reverses 
every 11 years around the maximum period of solar activity. By analyz- 
ing the GCR density measured with the anticoincidence guard on board 
ACE and Helios satellites and by ground-based neutron monitors (NMs), 
Richardson et al. [ 3 ] and Richardson [ 4 ] investigated GCR modulation by 
CIRs, which are formed at the leading edges of corotating high-speed solar 
wind streams originating in coronal holes on the solar surface. From statis- 
tical analyses of variations in many CIRs with and without the HCS, the 
above authors concluded that the IMF sector boundary (i.e., HCS) does not 
organize the GCR density, contrary to expectations from the simple drift 
model. On the other hand, they also found observational evidence that the 
amplitude of the CIR-related modulations of the GCR density around the 
solar minimum is ~ 50% larger in A>0 epochs than in A<0 epochs. The 
GCR modulation in CIRs is still an open question. 

These analyses were based on a one dimensional spatial distribution of 
the GCR density measured as the temporal variation of the GCR count rate 
recorded by each detector as it sweeps across the CIR. A single detector 
measurement of the density cannot, however, observe the spatial distribu- 
tion separately from the temporal variation of the density at a particular 
location in space. The spatial distribution sampled by a single detector also 
does not reflect the global distribution surrounding the detector, above and 
below the detector's viewing path. It is also affected by local irregularities 
in the solar wind. This can be serious, particularly when analyzing satellite 
measurements of low-energy (sub-GeV) GCRs. 

The spatial distribution of GCRs also can be inferred from the GCR 
spatial gradient in three dimensions, which provides us with additional 
information on the distribution around the detectors viewing paths. The 
diffusive flux equation describes the cosmic ray anisotropy in terms of the 
solar wind convection and the GCR density gradient. By inversely solving 
the diffusive flux equation, including the drift, Okazaki et al. [ 5 ] (hereafter 
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referred to as paper I) deduced the gradient from the anisotropy that is 
derived from the observation made by the Global Muon Detector Network 
(GMDN). The GMDN has a major response to 50GeV GCRs whose Larmor 
radii are as large as 0.2AU in 5nT IMF. They analyzed the three geocentric 
solar eclipic coordinate (GSE) components of the gradient observed in as- 
sociation with CIRs in a single solar rotation period (Carrington Rotation 
CR 2043) in 2006, around the solar activity minimum of A<0 epoch. It was 
found that the temporal evolution of the gradient before and after the HCS 
crossing is consistent with the prediction of the drift model, that the GCR 
distribution has a local maximum on the HCS [ : ][ 2 ]. As far as the gradi- 
ent in this particular rotation period concerned, on the other hand, there 
is no clear signature of the significant modulation in association with en- 
hancements in the solar wind velocity and the IMF magnitude in the CIR. 
Examining the drift effect solely from the gradient in a single rotation in 
A<0 epoch, however, is not easy, as the local maximum of the GCR density 
on the HCS is expected not only from the drift effect, but also from the 
minimum solar wind velocity on the HCS in a CIR. In the present paper, 
therefore, we examine the drift effect by analyzing the average gradient of 
multiple CIRs observed during the minimum periods of the solar activity 
in A<0 and A>0 epochs. 

2. Analysis and result 

2.1. Superposition analysis of CIRs at the HCS crossing 

We obtain the average feature of the GCR density gradient by superposing 
multiple CIRs at the timing of the HCS crossing in both A<0 and A>0 
epochs. Throughout this paper, we use the hourly mean solar wind velocity 
and IMF data from the Omni2 data for A>0 epoch and from ^4CE-Level2 
data for A<0 epoch, respectively 8,13 . We use ACE-Leve\2 data lagged by one 
hour, as a rough correction for the solar wind transit time between ACE and 
the Earth. By using IMF components (B x ,B y ,B z ) in GSE coordinates, the 
IMF sector is designated toward (T) if B x > B y and away (A) if B y > B x . 
We define the HCS crossing at the timing when the sector polarity (toward 
or away) designated by 23-hours central moving average of B x and B y 
alternates. We then choose HCS crossings associated by the enhancement of 
the solar wind velocity with the maximum velocity exceeding 500 km/s and 
the IMF magnitude with the maximum exceeding 10 nT. According to this 



a Omni2 data are available at http://omniweb.gsfc.nasa.gov/html/ow_data.html 
b ACE-Level2 data are available at http://www.srl.caltech.edu/ACE/ASC/ 
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Fig. 1. Sample HCSs superposed in this paper for A<0 (left) and A>0 (right) epochs. 
Each panel displays the solar wind velocity in a solar rotation period as a function of 
the universal time (day) on the horizontal axis. The time of the Earth's HCS crossing is 
indicated by the vertical solid line, while the away and toward IMF sectors are indicated 
by black and gray data points, respectively. Attached on the right side in each figure 
is the source surface synoptic charts (Radial 250) in the corresponding rotation by the 
Wilcox Solar Obscrvatory c 

selection criteria, we chose 42 CIRs observed in 2006-2008 for A<0 period 
and 38 CIRs in 1995-1998 for A>0 epoch, respectively. The number of CIRs 
superposed for A>0 period is limited mainly due to many data gaps in the 
Omni2 data. Figure 1 displays the solar wind velocity and the source surface 
synoptic charts (Radial 250 model) by the Wilcox Solar Observatory . The 
HCS crossings identified in this paper are also indicated by the vertical solid 
line in this figure. The average HCS tilt angle in A<0 epoch analyzed in 
this paper is 15.2°, while it is 9.0° in A>0 epoch. The superposed IMF and 
solar wind velocity are displayed in Figures 2a-e as functions of time for 
±5 days from the HCS crossing. Clear signatures of CIR around the timing 
of the HCS crossing, such as enhancements of the IMF magnitude and the 
solar wind velocity (Figures 2a and 2d) and reversals of the azimuthal solar 



c The source surface synoptic charts and the HCS tilt angle arc available at 
http: / / wso.stanford.edu/ 
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Fig. 2. From top to bottom, each panel shows the IMF magnitude (a), the GSE- 
longitude and latitude of the IMF (b-c), the solar wind velocity (d), the GSE- longitude 
of the solar wind (e) and three GSE components of the observed gradient (f-h). CIRs 
superposed at the HCS crossing in A<0 (left) and A>0 (right) epochs. Thin gray curves 
in each panel display individual CIRs superposed, while the black curve shows the aver- 
age as a function of time from the HCS crossing indicated by a vertical solid line. The 
dotted curves above and below the black curve bound the error of the average evaluated 
from the standard deviation of superposed CIRs. Only ±5 days period from the HCS 
crossing is plotted. The superposition at the HCS crossing from the northern (southern) 
hemisphere to the southern (northern) hemisphere is indicated by "N— >S" ("S— >N") at 
the top of each figure. 

wind flow angle (Figure 2e), are clearly seen in this figure. We also note 
that the superposed IMF longitude in away (toward) sector in Figure 2b 
is close to 135° (315°) and the average latitude in Figure 2c is close to 0° 
regardless the IMF sector. This implies that the nominal Parker spiral field 
with the spiral angle of 45° is a good approximation of the observed IMF 
as long as the average field is concerned. 

2.2. Analysis of GCR density gradient in three dimensions 

We deduce the GCR density gradient in three dimensions on hourly basis 
from the anisotropy observed with the network of muon detectors. Readers 
can refer to paper I for each muon detector in the network and also for the 
detail analysis method to derive the gradient from the observed anisotropy. 
In this paper, we analyze the average gradient in multiple CIRs observed in 
2006-2008 with the GMDN during the minimum period of the solar activity 
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in A<0 epoch. We also derive the average gradient in 1995-1998 in A>0 
epoch by using the observation with a couple of muon detectors at Nagoya 
and Hobart and examine the dependence of the gradient on the polarity of 
the solar magnetic field. 

We first correct the observed anisotropy vector in GSE coordinates for 
the solar wind convection anisotropy and the Compton-Getting anisotropy 
arising from the Earth's 30 km/s orbital motion. We then divide the cor- 
rected anisotropy vector £(t) into components parallel and perpendicular 
to the IMF as 

€(t) = €ii(t) + €-L(*)- (i) 

These components are given as 

t\\{t)=R L {t)a\\G™ ff (t), (2) 

Z±(t) = R L (t)[a ± G ± (t) - b(t) x Gi(i)], (3) 

where (t) and Gj_(t) are respectively the density gradient compo- 

nents parallel and perpendicular to the IMF, is the particle's effective 

Larmor radius of 50GeV protons in the observed IMF which is 0.2 AU in 
5nT IMF, h(t) is the unit vector pointing in the direction of the IMF, and 
ct|| and a± are the dimensionless mean free paths (A|| (t) and X±(t)) of GCR 
scattering by magnetic irregularities, defined as 

a„ =\\\/RL(t), (4) 

a± = X±/R L (t). (5) 

In this paper, we adopt an ad hoc choice of ctu = 7.2 and a± = 0.36 
which were also assumed in paper I. Qualitative features would remain 
similar for a wide range of reasonable parameters. Inversely solving Eqs.(2) 
and (3) for G™ ff (t) and Gj_(i), we obtain 

G ± (t) = Gf ff (t) + G^ H/t (t) (7) 

where, 

G? ff (t) = 7 — a 2 \ m t±(t) (8) 
(1 + a\)R L (t) 



December 24, 2009 13.32 



WSPC - Proceedings Trim Size: 9in x 6in ST-0083 



7 



and, 

G ?" , ''' ) = (l + aj W ) bX ^'>- < 9 ' 

Finally G(t) can be expressed in terms of two distinct contributions respec- 
tively from the diffusion and drift effects as, 

G(t) = {G° iff (t) + Gf ff (t)} + G^ ri/t (t), (10) 

where the first two terms in the bracket on the right hand side express the 
contribution from the parallel and perpendicular diffusions respectively, 
while the second term denotes the contribution from the drift effect. 




Fig. 3. The numerical solution of a simple drift model for 50 GeV GCRs. From top to 
bottom, each panel displays the GCR density and three GSE components of the gradient 
as a function of day in a solar rotation period in A<0 (left) and A>0 (right) epochs. 
The numerical values are shown by open circles. Contribution from the diffusion effect 
is shown by gray curve, while the contribution from the drift effect is displayed by black 
curve. The HCS crossing is indicated by a vertical black line. See paper I for the detail 
of this numerical model. 



Three GSE-components of the superposed G(t) are displayed in Figures 
2f-h, while the numerical prediction of a simple drift model are shown in 
Figure 3. This model assumes that the solar wind velocity and IMF magni- 
tude are constant without any enhancement seen as CIR signatures. This 
model is identical to that which paper I compared with the observation and 
readers can refer to paper I for more detail of this model. The observed gra- 
dient G(t) in Eq. (10) shows a large fluctuation due to the large fluctuation 
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of the IMF orientation (b x , b y , b z ). To show the systematic variation of each 
component, we filter this fluctuation with a central 23 hour moving average 
of hourly data. In Figure 2, the gradient in each CIR event is displayed 
by a gray curve, while the average gradient of all CIR events is plotted 
as the black curve. It is clear that the average GSE-z component of the 
gradient (G z ) in Figure 2h is negative (positive) above (below) the HCS in 
A<0 epoch, while it is positive (negative) above (below) the HCS in A>0 
epoch. This is consistent with the drift model prediction in Figure 3 at least 
qualitatively. We confirm this for the first time with the average gradient 
in multiple CIRs observed. It is also seen that the GSE-x component of 
the average gradient (G x ) remains negative in Figures 2f implying that the 
ambient radial density gradient is positive and the density increases away 
from the Sun in response to the outward convection by the solar wind. 
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Fig. 4. Contributions from the diffusion and drift effects to three GSE components of 
the average gradient G(t). From top to bottom, dotted curve in each panel displays the 
average G x , G y and G z in Figure 2, while gray and black curves show respectively the 



contributions to G(t) from the diffusion (G° iff (t) + G" lJJ (i)) and drift (G 



, Drift 



(*)) 



effects, in A<0 epoch (left) and A>0 epoch (right). G x and G y in this figure are obtained 
by averaging those in "N— >S" and "S— >N" in Figures 2f and 2g. We also calculate the 
average for G z in Figure 2h after reversing sign of G z in "S— »N" in Figure 2h. 



In Figure 2f, we also note the negative enhancement in G x observed 
following the HCS crossing. Such enhancements were first observed in both 
G x and G y in a sample rotation period in paper I and attributed to the 
drift effect causing the bidirectional density gradient toward the HCS in 
the ecliptic plane, together with the contribution from the ambient radial 
gradient away from the Sun. Similar enhancements are also seen in both 
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G x and G y in the drift model prediction in Figure 3. The observed gradi- 
ent in Figure 2, however, shows the enhancement only in G Xl but not in 
G y . This asymmetry is interpreted in terms of the contribution from the 
diffusion effect as follows. Figure 4 shows three components of the average 
G(t) by dotted curves, together with contributions from the diffusion ef- 
fect (G? iff (t) + G± iff (t)) by gray curves and the drift effect (G°" ft {t)) 

by black curves in Eq. (10). In A<0 epoch, G,, (t) due to the inward 

diffusion along the IMF is dominant as well as Gj_ *(t). G±" (t) shows 
the negative enhancement in both G x and G y , while G^ l ^(t) shows the 
negative enhancement in G Xl but the positive enhancement in G y . This is 
because of the geometry of the nominal Parker field in the GSE coordi- 
nate system. This asymmetric contribution from Gjj 3 (t) to G x and G y 
results in the asymmetry seen in Figure 2. This asymmetry, on the other 
hand, is less evident in A>0 epoch, because of another contribution from 
G® 1 ™ (t) due to the inward perpendicular diffusion which results in nega- 
tive enhancements in both G x and G y . These contributions to the observed 
G(t) from the parallel and perpendicular diffusions are not clarified by a 
sample rotation analysis in paper I and first revealed in this paper by ana- 
lyzing the average gradient in multiple CIR events. We are now preparing 
for a quantitative analyses of each contribution from G® 1 ^ (t), G 1 ^ 1 ^ (t) 

and Gf i/4 (i). 




Fig. 5. GCR density superposed at the HCS crossing in A<0 (left) and A>0 (right) 
epochs. This figure is plotted in the same format as Figure 2, except that the derrived 
density is normalized to % at the HCS crossing to clarify the variation associated with 
the HCS crossing. 



In Figure 5, we finally display the GCR density which is derived as 
the omnidirectional intensity, together with the anisotropy vector, from the 
best-fitting to the observed muon count rates in various directional channels 
(see paper I). In this figure, we normalize the density to % at the HCS 
crossing to clarify the variation associated with the HCS crossing. It is seen 
that the amplitude of the average variation is small and the fluctuation is 
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large even for 50 GeV GCRs recorded with muon detectors, easily masking 
the average systematic variation seen in figure 3 by the numerical model. 
This difficulty due to the fluctuation can be serious, particularly when an- 
alyzing low-energy GCRs, which have smaller Larmor radii and are more 
sensitive to small-scale structures in the solar wind. The variation in this 
figure is much less significant than the average gradient variation in Figure 

2. This implies how difficult is to deduce the GCR spatial distribution solely 
from single detector measurement of the GCR density[ 3 ][ 4 ]. 

3. Discussions 

In table 1, we compare the average GCR gradient and density with the 
model prediction. We first note that the magnitude of G z is 50 % larger 
in A>0 epoch in accordance with the model prediction. The magnitude of 
negative G x , on the other hand, is larger in A>0 epoch, while the model 
predicts G x almost twice in A<0 epoch in an oppsite sense. We also note 
the average density (Jo) being 50 % lower in A>0 epoch. As we normalized 
Jo to % at the HCS crossing, this implies the amplitude of the density 
variation before and after the HCS crossing is 50 % larger in A>0 epoch, in 
accordance with Richardson et al. [ 4 ] reporting that the amplitude of the 
recurrent GCR intensity variations in solar rotation periods are 50 % larger 
in A>0 epoch than in A<0 epoch. 

In concluding anything from these observed features about the drift 
model prediction, however, we should be aware of that the gradient in 
A>0 epoch presented in this paper might be biased, because it is derived 
from an incomplete network before the GMDN started operation in March, 
2006. To check such a possible bias, we analyzed A<0 epoch by using only 
two detectors (at Nagoya and Hobart) used for A>0 in the GMDN and 
confirmed that the derived gradient can be 10 ~ 30 % larger than that 
derived from the entire GMDN depending on values of a\\ and a± assumed. 
In order to precisely clarify the difference of the gradient in A<0 and A>0 
epochs, therefore, we need to improve the analysis method further and 
correct the observed gradient in A>0 for this bias effect. 
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Table 1. GCR density gradient (G x ,G y ,Gz), density (lo), solar wind velocity (Vsw) 
and IMF magnitude (B) averaged over a period in Figure 2 in A<0 and A>0 epochs. 
The upper lines for the gradient and density represent the observed average, while the 
lower lines represent the model prediction. The average G z is calculated after revers- 
ing the sign of G z in toward sectors. The period between ±1 day from the HCS cross- 
ing is excluded from the average calculation. Errors for the observed values are de- 
duced from the standard deviation of CIRs used for the calculation, while errors for the 
model prediction are derived from the standard deviation of data points in Figure 3. 
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